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We describe a model independent lattice QCD method for determining the deviation from unity for /iai(1), 
the B -^ D*t.u form factor at zero recoil. We extend the double ratio method previously used to determine the 
B -^ Dlv form factor. The bulk of statistical and systematic errors cancel in the double ratios we consider, 
yielding form factors which promise to reduce present theoretical uncertainties in the determination of |Vcf,|. We 
present results from a prototype calculation at a single lattice spacing corresponding to /3 = 5.7. 
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1. Introduction 

The form factor Ha^ (1) parameterizes hadronic 
matrix elements in B — > D*ii/ decays. Its theo- 
retical determination is necessary in order to ex- 
tract the CKM matrix element \Vcb\ from the ex- 
perimental decay ratepl, extrapolated to zero re- 
coil. 
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Heavy quark symmetry constrains this form 
factor[p|. Up to radiative corrections, it has 
deviations from unity beginning at order I/tiq 
in an expansion in inverse powers of the quark 
masses 0. It is exactly one in the infinite mass 
limit. We write: 



/i^,(l) = 77A 1 - (5i/„2 + 0(1/ 



m|) 



(2) 



where the heavy quark expansion in Heavy Quark 
Effective Theory (HQET) is within brackets, and 
rjA denotes radiative corrections in the HQET-to- 
QCD matchingH). 

Previous calculations of 5-, /,„2 have relied on 
quark models or sum rules estimates [y. We de- 
termine (5x/„j2 ; in principle model independently, 
using lattice QCD. Our method extends our pre- 
vious work in determining the zero-recoil form 
factor in B —>^ Diu decays p. 

* Presenter. 
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. Procedure 










Consider double ratios, 






E 


jfr^'w - 










\B) {B\J, 
\D) {B\J^ 


D) 




T/2»t»0 


\B) 



(3) 



of lattice three-point functions, C^'^'^^{t) ~ 
\Xd{T/2) Jp,{t) Xsi^))- E)ouble ratios are con- 
structed to be identically one when the "charm" 
and "bottom" quarks are of equal mass. The bulk 
of statistical and systematic uncertainties cancel 
in such ratios 1^. 

We need three double ratios: 
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The lattice-to-HQET matching coefficients, p? , 
are near unity for typical lattice spacings0. 
Right-hand expressions in the equations above 
are to be interpreted within HQET. Known nor- 
malizations for elastic vector-current matrix ele- 
ments were used to simplify the first two ratios. 
These normalizations are obtained nonperturba- 
tively in our numerical work. 
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Table 1 

Ratios and their l/mg coefficients. 



Ratio 


coefficient c- ' 
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All three ratios have quark mass dependence 
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where, A„q = (^-^ - ^j . Table |i| displays 

(2) 

our notation for the q \ Note that the B -^ 
D*£i> coefficient contains a linear combination of 
the other two coefHcients. The form of this coef- 
ficient is derived by substituting the mass depen- 
dence ||], 

-A (^ ^''^ ^ 



-51/ 



(8) 



27Ti(,y Arricmi, 

into the expression for R^^^ shown in Eq. ^ 
Our procedure for determining ^.^^^(l) in _B ^ 

(2'] 

D*(,v decays is: Extract c\ by studying the mass 
dependence of the double ratios. Solve for ip^ ly 
and A. Substitute these values and the values we 
determine for rric and mf, into the expression for 
5ii^2 shown in Eq. g. We then match to QCD 

and determine /i^i(l) as in Eq. ||. 

3. Prototype Calculation at /3 = 5.7 

Many of the numerical details in this study are 
common to our study of B -^ Div matrix ele- 
ments [|6|. We note in particular: 

• We use a subset of 200 P — 5.7 quenched gauge 
configurations on a 12'^ x 24 lattice. 

• We use the Sheikholeslami-Wohlert quark ac- 
tion with a tadpole-improved tree-level coeffi- 
cient, csw ~ 1.57. Results arc interpreted within 
the Fcrmilab heavy quark formalism[g|. 

• We study double ratios for nine combinations 
of heavy quarks with bare masses corresponding 
to Kh e {0.125,0.119,0.110,0.100,0.089,0.062}. 



Quark masses range around both the charm and 
bottom masses. 

• Statistical errors were obtained using a single- 
elimination jackknife procedure. 

• The physical tree-level charm and bottom 
quark masses, rric and TTib, were determined by 
adjusting bare mass inputs and demanding that 
calculated meson kinetic masses match the phys- 
ical D and B meson masses. 

• Matching factors pj are only known to one 
loop order. For consistency, r]A is truncated to 
one loop order. Matching factors are computed 
using the V scheme coupling. We use BLM 
matching scales which account for /^o^s contri- 
butions. 

Results in this paper have the spectator quark 
mass fixed near the strange quark mass. In our 
B — > Dtv study, we checked the dependence upon 
the spectator mass for Ry^^. Values in the chi- 
ral limit were consistent with those for the strange 
quark, but with statistical errors which were twice 
as large. We anticipate similar chiral behavior for 
the other two ratios we use in this study. Hence, 
we expect insignificant differences in c- in the 
chiral limit, and similar increases in statistical er- 
rors. Here, we account for the uncertainty of not 
performing the chiral extrapolation by doubling 
our statistical errors. 

Figure IT] shows the heavy quark mass depen- 
dence we find for \ R^^^' ■ The quality of these 

results are representative of our results for the 
other two ratios. The two points to the left in 
the figure have large statistical errors. These de- 
cays involve the heaviest quark masses and suffer 
from well-known signal-to-noise problems. A fit 
to the functional form given in Eq. is shown in 
the figure. Its y-intercept shows this c^^' deter- 
mination has greater than 4cr significance. 

We find this and the other two coefficients are 
of 0{Kqcd), as expected. The values we ob- 
tain are broadly consistent with previous esti- 
mates [^. Quantitative comparisons may be mis- 
leading, however, since uncertainties in previous 
estimates are difficult to ascertain. 
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Figure 1. Calculated mass dependence for 

Jr^-^°' at /3 = 5.7. The solid line is a fit 

to the expected mass dependence given in Eq. H. 
Dashed lines show 1-a statistical errors from the 
fit. The burst indicates the fitted ratio at the 
physical combination of quark masses. 



the decay constants fo and /b and oi B —^ ■nf.u 
matrix elements we find lattice artifacts are un- 
der control [g 10 1 . We anticipate no large lattice 
artifacts in this study. We adjust the quark ac- 
tion and currents to match Vttt-o terms of the 
QCD heavy quark expansion[p|. Contributions 
to (Si /,„2 arise solely from these terms in the dou- 
ble ratio methodW. We note that any remaining 
cutoff dependence may be removed by repeating 
the calculation for additional lattice spacings and 
taking the continuum limit. 

• Uncertainty due to the quenched approxima- 
tion. A full quantitative estimate of the error due 
to quenching must await an unqucnched deter- 
mination of ft-Ai- Note, however, the uncertainty 
due to quenching affects the deviation of the form 
factor from unity. The quenching uncertainty in 
this deviation is commonly believed to be perhaps 
20%. 



4. Determination of hAi(l) 

We present a preliminary determination of 
/lAi (1) using our prototype (3 = 5.7 study to illus- 
trate the precision we expect in a complete study: 

/lAi (1) = 0.935 ± 0.022( tlf^l) ± 0.008 ± 0.020. 

Sources of uncertainties are, respectively: 

• Statistical. We double statistical errors to ac- 
count for not having extrapolated to the down 
quark mass. We must still check the chiral be- 
havior of all three double ratios used in the de- 
termination of Hai ■ 

• Tuning of quark masses rric and rub. We es- 
timate 10% and 13% uncertainties in our charm 
and bottom masses. 

• Unknown radiative corrections beyond 1-loop. 
We estimate this uncertainty by varying the 1- 
loop coefficients by 20%. 

• Undetermined 0{}/m%) corrections to /lyi^(l). 

We use the relative sizes of the ^Im?^ terms de- 
termined in our fits to estimate neglected power 
corrections. 

Two important sources of systematic uncertainty 
remain to be evaluated fully: 

• Lattice spacing dependence. In our studies of 
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